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Abstract 

We consider theories which explain the flatness of the power spectrum of scalar 
perturbations in the Universe by conformal invariance, such as conformal rolling model 
and Galilean Genesis. We show that to the leading non-linear order, perturbations in 
all models from this class behave in one and the same way, at least if the energy 
density of the relevant fields is small compared to the total energy density (spectator 
approximation). We then turn to the intrinsic non-Gaussianities in these models (as 
opposed to non-Gaussianities that may be generated during subsequent evolution). The 
intrinsic bispectrum vanishes, so we perform the complete calculation of the trispectrum 
and compare it with the trispecta of local forms in various limits. The most peculiar 
feature of our trispectrum is a (fairly mild) singularity in the limit where two momenta 
are equal in absolute value and opposite in direction (folded limit). Generically, the 
intrinsic non-Gaussianity can be of detectable size. 



1 Introduction 



Scalar perturbations in the Universe may well originate from inflation In that case, 
approximate flatness of their power spectrum is due to the approximate de Sitter symmetry 
of inflating background. As pointed out in Ref. 0, a possible alternative to the de Sitter 
symmetry in this context is conformal invariance. Concrete models employing conformal 
invariance, instead of the de Sitter symmetry, for explaining the flat scalar spectrum include 
conformal rolling and Galilean Genesis [Q. Despite substantially different motivations 
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and dynamical features, the latter models end up in one and the same mechanism that 
generates scalar perturbations. From this prospective, the theory boils down to a model 
with two scalar fields p and and the Lagrangian 

L = L, + {d,ef , (1) 

where Lp governs the dynamics of p. Under scaling transformations, the field p scales as 
p{x) — )■ Xp{Xx)] it provides for a non-trivial background pc- The field scales as 0(a;) — )■ 
0(Aa;); the perturbations of the field serve as predecessors of the adiabatic perturbations. 
One makes sure that at the time the perturbations of are generated, gravitational effects 
on the dynamics of the two fields are irrelevant (this is achieved in different ways in Ref. 
and Ref. [^) and assumes that the background field pc is spatially homogeneous. Then 
conformal invariance of the field equation for p implies that 

Pc{xq) = — - , xo < , (2) 

Xo 

where Xq = 1] is conformal time in the conformal rolling model and a;o = ^ is cosmic time 
in the Galilean Genesis scenario. In such a background, the field behaves exactly in the 
same way as massless scalar field (minimally coupled to gravity) in the de Sitter space-time. 
The modes of its perturbations about a certain background value 0, 

e{x) = 0(x) - , 

start off in the WKB regime and freeze out when k\xo\ ~ 1, where k is (conformal) mo- 
mentum. Assuming that the field 6{x) is originally in its vacuum state, one obtains, at 
the linearized level, the fiat power spectrum of the Gaussian field ^(k) at late timesQ, when 
k\xo\ <^ 1. The conformal stage ends up at some point, and the 6'-perturbations are repro- 
cessed into adiabatic ones at much later epoch via, e.g., curvaton P or modulated decay ^ 
mechanism. 

A common feature of the conformal models of Refs. Q is the existence of the per- 
turbations of the field p about the background (|]). The Lagrangian Lp contains a small 
parameter, call it h, then 6p (x h [h"^ is quartic self-coupling in the conformal rolling model; 
we relate the parameter h to the parameters of the galileon Lagrangian in the Galilean Gen- 
esis scenario in Section pTTI , see Eq. (H)). As discussed in Refs. ^, the perturbations 6p 



have rather peculiar properties. Nevertheless, as we review in Section |2.1| , these perturba- 
tions are exactly the same (modulo overall amplitude) in the two, apparently very different 
models. Furthermore, in Section |2.2| we give a general argument showing that the properties 



"'^ Deviation from exact conformal invariance naturally gives rise to the tilt in the power spectrum, see, 
e.g., Ref. i- 
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of the linearized perturbations 6p are common to a large class of models employing conformal 
invar iance. 

There is a subtlety here. Strictly speaking, our observation is valid in a theory without 
gravity or in the case when the energy density of the field p is small compared to the total 
energy density in the Universe (spectator approximation). Though is is likely that our results 
are valid in much more general setting, the corresponding analysis is yet to be done. We 
proceed in this paper by neglecting gravity altogether. 

So, it is of interest to study the effects of the perturbations 5p on the field 6 and, in 
the end, on the adiabatic perturbations. There are at least two of these effects, namely, 
the statistical anisotropy and non-Gaussianity. Of course, the non-Gaussianity may well be 
generated also at the time the ^-perturbations are converted into the adiabatic ones. In 
this respect the conformal models are nothing special as compared to inflationary theories 
equipped with the curvaton or modulated decay mechanism; for this reason we are going to 
disregard this conversion-related non-Gaussianity. What we are interested in is the intrinsic 
non-Gaussianity, which is due to the interaction of the field 9 with perturbations 5p. 

The resulting phenomenology depends strongly on what happens to the field 9 after the 
conformal stage @ ends up. One option is that the field 6'(k) starts evolving again, and its 
evolution proceeds until it becomes superhorizon in the conventional sense. This option is 
fairly natural in the conformal rolling model of Ref. and more contrived in the Galilean 
Genesi^ of Ref. [Q. Both statistical anisotropy and non-Gaussianity generated in this case 
are studied in Ref. p. 

Here we consider the opposite case, i.e., we assume that the field 9 does not evolve after 
the end of the conformal stage. This option is particularly natural in the Galilean Genesis, 
but it is not contrived in the conformal rolling scenario either. In models from this sub-class, 
the adiabatic perturbations inherit the statistical properties of ^-perturbations that exist 
already at the late conformal stage, when fc|a;o| ^ 1. The statistical anisotropy in the field 
9 is induced by long-ranged perturbations 5 p. This effect has been studied in Ref. with 
the result that the leading statistical anisotropy in the power spectrum of 6', and hence of 
adiabatic perturbation has the quadrupole form. 



where Vo{k) is isotropic (and nearly fiat) power spectrum obtained at the linearized level, Ui 

^In the context of the Gahlean Genesis, an intermediate stage of the evolution of 6 may occur provided 
that the effective scale factor p^ir) is a non-trivial function of the galileon field tt, such that p cx e'^ at tt 
smaller than some value ttq and p = const at tt > ttq. Then the field 9 feels the background at early times, 
when TT < TTo, and temporarily gets frozen out when k\xo\ ^ 1, as discussed in the text. If the opposite regime 
TT > ttq sets in later, but still at some sufficiently early time when the space-time is nearly Minkowskian, 
the field 0(k) indeed starts to oscillate again at that time. These oscillations terminate when the Hubble 
parameter becomes large enough and the field 9{k) exits the horizon. 
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and Wij are unit 3- vector and unit traceless 3-tensor of a general form, Hq is the present value 
of the Hubble parameter, the parameter ci is of order 1 and C2 is positive and logarithmically 
enhanced. 

The main purpose of this paper is to study the intrinsic (as opposed to conversion-related) 
non-Gaussianity in this sub-class of models. In the absence of the cubic self-interaction of the 
field 6, the intrinsic bispectrum vanishes, so we have to consider the intrinsic trispectrum. 
Perhaps the most striking feature of the trispectrum is the singularity in the limit where 
two momenta are equal in absolute value and have opposite directions (folded limit, in 
nomenclature of Ref. [0). The singular part of the connected four-point function has been 
calculated in our earlier paper [O with the result 



(CkiCk2Ck3Ck4) = const -5 i^k. 



, j=i 

ki2 = ki + k2 

-1 



ki2ki 

ki2ki 



ki2k3 

ki2k2, 



(3) 



i.e., the trispectrum blows up as kT^2 ■ This is in contrast to trispectra obtained in single- 



field inflationary models |T0, O, O, Ol |T5[, and, indeed, there are general arguments |T3 



showing that in these models, the four-point function is finite in the limit ki2 — ?• 0. The 
singularity in the four-point function (^ is due to the enhancement of the perturbations 5p 
at low momenta, see the discussion of the infrared properties of the conformal models in 
Ref. [^. We will see in Section 173 that the most relevant features of the trispectrum are 
captured by its singular part. 

Whether the intrinsic non-Gaussianity dominates over the conversion-related one, and 
whether the former is detectable depends on both the underlying conformal model and the 
mechanism that converts the 6'-perturbations into the adiabatic ones. Generically, the linear 
order relationship between C, and 6 is ( = r6/Q, where r < 1 is the dilution factor. With 
our normalization (|l]), (Q), the power spectrum of ^-perturbations is Ve = l/(47r^). The 
background value G cannot be estimated in a model independent way. It can be as large 
as ~ 10^, so that the correct adiabatic amplitude is obtained at r ~ 1. This can be 
the case, e.g., in the Galilean Genesis model. In such a situation, the conversion- related 
non-Gaussianities induced by the curvation or modulated decay mechanism are fairly small 
{Inl, Qnl are roughly of order 1 [|16|, [T^, 0]). On the other hand, the size of the intrinsic 
non-Gaussianity t^L (see Eq. ( ^Bf ) for its definition) is governed by the amplitude of the 
perturbations 5p, so that 

see Eq. (^) for numerical coefficient. There is no general reason to expect that the parameter 
h is particularly small, except for the mild requirement /i <^ 1 ensuring the self-consistency of 
the conformal scenario. This shows that the intrinsic non-Gaussianity may well be dominant 



4 



and detectable (say, t^L ^ 10^), though it size cannot be predicted because of our ignorance 
of the value of the parameter h. The situation is more subtle in the conformal rolling model; 
we consider this point in Section with the result that dominant and detectable intrinsic 
non-Gaussianity is possible, but not generic. 

This paper is organized as follows. To set the stage, we consider in Section |2.1| the lin- 
earized perturbations 6 and 6p in the Galilean Genesis scenario and compare the resulting 
expressions with those obtained in the conformal rolling model. This Section contains noth- 
ing new as compared to Refs. 0, ^ ; the main point is to show that the perturbations 6 and 
6p are identically the same in the two scenarios. In Section p.2| we give a general argument 
showing that the properties of 6p are uniquely determined by conformal invariance (modulo 
the overall constant amplitude). Hence, the peculiarities of the adiabatic perturbations that 
we study in this paper, as well as the results of Refs. |1T[, are common to the whole class 
of conformal mechanisms. We turn to the non-Gaussianity in Section ^ where we perform 
the complete calculation of the intrinsic trispectrum. We confirm our earlier result (|^) con- 
cerning the singular part of the trispectrum. We then consider various limits, following the 
nomenclature of Refs. [0, and compare them with local models [^, We conclude 
in Section |. Some details of our calculations are collected in Appendix. 



2 Perturbations in conformal scenarios 
2.1 Galilean Genesis vs. conformal rolling 



The galileon model has been introduced in Ref. |2^. The rolling galileon serves as the 
field p entering (|^) and (0). In the Minkowski space-time, the Lagrangian of the simplest 
conformally-invariant version of the model is (mostly negative signature) 



= -pe'^d^nd^n + ^d^ndWn + ^(^^ttS^tt)^ , (4) 



where □ = df^d^. The field equation in Minkowski space-time admits the homogeneous 
solution TTc such that 

where xo = t and 
Hence, one defines 

p = i/*e" (6) 
so that the background solution is given by @. 
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The quadratic action for perturbations about this solution is 

^* J 1^0 ^0 J 

It is convenient to introduce the variable Sp = —Xq^Stt. It follows from @ that 5p is the 
perturbation of the field p about the background pc = —Xq^. Its action reads 



S, 



Sp 



2P 



d X 



The field equation has a simple form, 



{d,5pf-{d,5pf + \{5p) 

Xn 
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d'^dp - A6p - —5p = 



Xn 



(7) 



Its properly normalized solution for given momentum k is 



6p = h ■ e 



ikx 




4^ V i-^^o) ■ + h.c. 



(8) 



where bI. and 5^ are creation and annihilation operators obeying the standard commuta- 
tional relation [B^, = 5(k — q), and 



h 



V2f 



(9) 



The theory is weakly coupled in the relevant range of momenta provided that h '^1. 

The mode (|^) oscillates at early times, when fc|xo| ^ 1, while at late times it behaves as 
follows. 



Sp = h ■ e 



ikx 



1 



4^3/2 ^5/23,2 



5k + h.c. 



This behavior can be interpreted as local time shift Indeed, for time-shifted galileon 

solution we have 



e'^=(l + (57r) 



1 



1 6x0 

+ 



H^,{xo + 6xo) H^xq H^xl 



2 • 



Hence, 



5xq = —xoStt = xlSp . 
Thus, (5xo(x) is independent of time at late times, and its power spectrum is red, 

9/^2 1 



47r2 /c2 



(10) 
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In the Galilean Genesis scenario, the field G is introduced as an additional field precisely 
for the purpose of generating the scalar perturbations. By conformal invariance, its quadratic 
Lagrangian has the form of the second term in (|I|). In the background pc = —Xq^, its modes 
are 

e = e'kx . __^(_a;Q)3/2^{i) _ ^ ^ ^^^^ 



TT 



where Al^ and is another set of creation and annihilation operators. At early times, the 
field 6 is in the WKB regime, while for k\xo\ ^ 1 the mode stays constant in time. The 
resulting late-time power spectrum is fiat. 

This result is valid at the linearized level. The lowest order interaction between 6 and 6p is 
described by the interaction Hamiltonian, whose density is (in the interaction picture) 

Hi = -Lint = -Pc Sp {d^ef . (13) 

It is this interaction that is responsible for the intrinsic non-Gaussianity which we study in 
Section ||. 

In the above discussion we neglected gravity effects. This is legitimate in the Galilean 
Genesis scenario at the early stage of Genesis, when the energy density of the galileon field 
is small, while the energy density of other fields is assumed to vanish. The latter assumption 
is relaxed in the conformal rolling scenario whose main ingredient is a complex scalar field 
conformally coupled to gravity. Unlike in the galileon case, the scalar potential does not 
vanish; it is assumed to be negative and is quartic by conformal symmetry, V{(f)) = — 
One makes use of the parametrization 

^ = ^xp (M) . (14) 



Then the background field pc rolls down the potential according to (^, where Xq = r] is now 
conformal time [Q. The Lagrangian for G coincides with the second term in (|Tp, while the 
perturbations 6p are again governed by Eq. (^. The properly normalized solutions for the 
field 6 and perturbations 6p coincide with ( [Til) and (^, respectively, with h being now the 
quartic self-coupling and Xq = rj. So, the dynamics of perturbations in the conformal rolling 
scenario is identical to that in the Galilean Genesis. 



2.2 General argument 

Let us see that the form of Eq. (|^ that governs the perturbations Sp at the linearized level is 
completely determined by conformal invariance. Namely, let us consider any weakly coupled 



7 



theory, in which the classical field equation for p is second order in derivatives, the action for 
p is local, invariant under space-time translations and spatial rotations and invariant under 
scaling 

p{x) — )■ Xp{\x) 

and inversion 

One designs a theory in such a way that it admits the runaway solution pc = —Xq^ (one 
can always set the overall constant in p^ equal to 1 by field redefinition). This solution is 
invariant under both scaling and inversion. So, the quadratic action for perturbations 6p 
about this solution must also be invariant. Let us write for the quadratic action 

S^^^ = h-^ j d^x 6pC6p, (15) 

where C is second order differential operator, and h is some constant whose choice is specified 
below. Since the background depends only on time, the operator C does not contain spatial 
coordinates explicitly, but may contain time. Also, it is invariant under spatial rotations. 

Invariance under scaling implies that modulo overall constant, the part of C that involves 
derivatives is 

CD-d'o+v^A. (16) 

where Vs is independent of time. We can choose the constant h in (|T5|) in such a way that 
the term with two time derivatives enters with the coefficient —1. Scale invariance alone is 
insufficient^ for obtaining Vg = 1. However, the requirement of invariance under inversion 
uniquely specifies Vg = 1- The term without derivatives is almost uniquely determined from 
the requirement of invariance under scaling and inversion, 

£ D ^ (17) 

where c is yet undetermined constant. To complete the argument we notice that the in- 
variance of the original action under time translations implies that Sp = pc = Xg ^ must be 
a solution to equation C6p = 0. This gives c = 6. Thus, the whole quadratic action for 
perturbations is uniquely determined by conformal invariance, modulo an overall constant 
factor, and the linearized equation for 6p has one and the same form in the whole class of 
conformally-invariant models with the Lagrangians of the general form (|l]). To the leading 
non-linear order, the properties of the 6'-perturbations are identical in these models, as they 
are governed by the interaction Hamiltonian (|13[). 

■^A straightforward way to see this is to consider a modification of the gahleon theory in which the second 
and third terms in (H) enter with unrelated coefficients. 
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3 Trispectrum 



3.1 Generalities 



In models with the Lagrangians of the form ([^), the bispectrum of the field 9 vanishes. To 
calculate the trispectrum we make use of the in- in formalism, cf. Ref. [Q. We use the 
shorthand notation 

ei ^ e{x)e{y)e{z)e{w) , 

with understanding that we are interested in the formal limit 

Xq = yo = zo = Wq . (18) 
Then the four-point function reads 



Texp I i / dxoHj 



Texp 



cLxqHi 



where subscript / refers to interaction picture and the interaction Hamiltonian density is 
given by (|T^. Of course, we are going to calculate the connected part of the four-point 
function. 

To this end we need the two-point functions of the linear fields 9(i) and 5p(/), 

= D{x, y) , {6p(i){x)6p(i){y)) = —Dp{x, y) . 



By making use of (0) and (||) we find 

1 



D{x,y) 
Dp{x,y) 



327r2 
1 



d^k {x,y,f/^ H^^]li-kxo)H^%{-kyo)e'''^''-^^ 
d^k {xoyofl^ HfJ^{-kxo)Hfl>^{-ky,y''^^-^^ 



(19) 
(20) 



32^2 / - - V— ^ ^^5/2 

Both of these pairing functions satisfy 

D{p){,x,y) = Dl){y,x) . 

We also need the (anti)T-product of the field 6 p. We write 
— 2 — 

Tccy + T^y = j^[{T{5p(i){x)5p(i){y))) + {T{5p(i){x)5p(i){y)))] = 

= e(xo - yo)Df,{y, x) + Q{yo - xo)Dp{x, y) + Q{xq - yo)Dp{x, y) + 0(?/o - xo)Dp{y, x) 
= Dpix, y) + Dp{y, x) = 2Re[D,(x, y)] (21) 
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and 



Txy - T^y = -e{xo - yo) y) - Dp{y, x)] = -2ie{xo - yo)lm[Dp{x, y)] . (22) 

The connected four-point function is a sum of the three terms, one of which has the following 
form, 



i I d^x'Hiix) I efj) I -i I d^x"Ui{x"^ 



2h' 



d^x'd^x" 



XqXq 



9' x)d' D{x', y)&lD{z, x")&lD{w, x")Dp{x', x'' 



X z 
y ^ w 

and the two others are 



+ {y^z} + {y^w} 



(23) 



- / d^x'd''x"THi{x')Hi{x") I 9\,^ ) + { ^f/) , 2 



d^x'd*x"THi{x')Hi{x") 



d oc d 

XqXq 



d'^D{x\ x)d'^D{x\ y)d'lD{x", z)d'lD{x", w)T,,,„ + 



+ d'D{x, x')d'D{y, x')d'lD{z, x")d'lD{w, x")T,'." 



+ {y^z} + {y^w}. (24) 



In the 3-dimensional momentum representation, these expressions reduce to the integrals 
over dxQ and dx^. Our definition of the n-point function in the momentum representation is 

ill ^k.) = i27if6 5^ k, G„(ki, . . . , k„) = / n {d'x.e-''^'^' } G„(xi, . . . , x„), 

i=l \i=l J i=l 



which corresponds to 



e(k) = j , 



d^xe-'^'^e^x) . 



With our convention, the power spectrum is related to the two-point function by 

(e(k)e(kO) = (2vr)^5(k + k')^P,(fc) 



and at the linearized level Ve is given by (|T2| 
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3.2 Leading singularity 



Let us calculate the contribution to the trispectrum due to the term explicitly written in 
square brackets in (^^. In the limit (|18D it involves the following combination 



d'Dix',^)d'Dix\y) 



S'kiS'k 



^giki(x'-x)+ik2(x'-y) 



1 2 /"rf'M^ ,:.,rx'-xw.k,rx'^va^^^^^^^^^^^^^)^^g^^^^^^^^^^^^) ^^5) 



(327r2)2 7r2 



1^11^2 



2giki(x'-x)+ik2(x'-y) 



where 



R{h, k2, fci2, = cos(e(fci + k2)) ■ [ekMkf^ - {ki + k^f) + kl + kl - kl^] 



+ sin(^(/ci + fca)) ■ ^(A;i + k2){kl + kj - k 



12) 



(26) 



and 



Q(fci, k2^ k,2, = sin(e(fci + A;2)) ■ [ekMkl, - {k, + k^f) + kl + kl - kl^] 

- cos{i{ki + k2))-i{ki + k2){kl + kl-kl2) . (27) 

with ki2 = kl + k2. So, we obtain for the contribution under study (in a certain sense this 
is the leading contribution, hence the notation) 



G'4a.c.)(ki,k2,k3,k4) 



oo j,3i,3l,3i,3 



X 



Q{ki, k2, ki2, OQih, h, ki2, x) [Y5/2ih20Y5/2ih2X) + J5/2{ki2i)J-,/2{ki2X 



where 



X 



and we made use of the fact that ks + k4 = — ki2. To regularize integrals at infinity we 
insert exp(— £^) exp(— £%) into the integrand and take the limit e — )■ +0 in the end of the 
calculation. This is equivalent to the replacement Xq — )■ Xo(l — is) which is the standard 
prescription for calculating vacuum expectation values in the interaction picture (see, e.g.. 
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Ref. p3|). At ^ — > all integrals are converging. We observe that the integrals over ^ and 
X factor out and obtain 



G4(i.c.)(ki,k2,k3,k4) 



[3^(A;i, ^2, fci2)3^(/i;3, ^4, A;i2) + J{ki, k2, h2)Jih, h, ki2)] 



where 



y{ki,k2,ki2) 



and 



-y=Q{ki, /C2, ki2, O^^MkuO 



2"^/^^^ ["^(^1 ^2)^ "^{kl + kl)kl2 ^^12] 



(2^ 



(29) 



k2, k 



12 



-J=Q{ki, /C2, /Cl2, 0'>^5/2(A;i20 



A;i2(A;L-3(A;i-A;2)')(A;i + A;2) 



/2^A;;^/' 



+ (3(A;^ - klf - 2{kl + kl)kl^ - k^^) arctanh 



ki2 



ki + A;2 



(30) 



The expression (^Uj) is singular in the folded limit ki2 — )■ 0, while J^{ki, k2, k^) vanishes in 
this limit. The crossing terms in (|23|) , as well as the terms (p^ are regular as ki2 — ?■ (the 
latter property is established by direct calculation in Appendix). It is convenient to define 
the singular part of 3^ in the following way, 

ysikl, k2, k,2) = - - 2(^1' + ^2)^L + k{2] , 



■^12 

SO that 

y = ys + yr. yr = 

Then the singular part of the four-point function is 

Tlh^ 1 



^r.3/2 

2 12 



G4(^)(ki,k2,k3,k4) 



09 j,3l3l3l3 



3^3(^1,^2,^12)3^5(^3,^4,^12) + 2 permutations (31) 



In the limit ki2 — )■ we recover the result (^, which can be written in a symmetric form, 



G4(,)(ki,k2,k3,k4) 



1 



|fcl2-i'0 



16 ki2Q^P* 



1 - 3 



ki2Q 

A:i2<5 



1-3 



ki2P 



A:i2P 
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where 

^ _ ki - ka ^ _ kg - k4 

^" 2 ' ^"2 
Note that the terms of order 1 / fc^g ^md 1 / , which one could off hand expect from (|29|) , 
cancel out. Thus, the field 6 has rather mild infrared behavior, even though it interacts with 
infrared-enhanced modes of 6p (the power spectrum of 6p is red, see ([T0|) ). The reason for 
this property is discussed in Ref. 0. 



3.3 Shapes 



The calculation of the contribution ( plD is not so straightforward. We perform this calcu- 
lation in Appendix, where we also present the complete result for the trispectrum. The 
definition of the trispectrum T in our case is 



G'4(ki,k2,k3,k4) 



i=l 



T{ki, k2, ks, ki, ki2, k 



14; 



where ki4 = ki + k4. The combinations ki^ 
independent, since 

kl3 = ^24 = 



|ki + kg I and k 



24 



k4| are not 



k-^ ~\~ ~\~ ~\~ k 



^14 



It is convenient for the purpose of illustration to decompose the trispectrum into the part 
Ts, which is singular as either ki2 — )■ or /cia — )■ or /ci4 — )■ 0, and the regular part %, 

r = Ts + Tr. 



The singular part is obtained^ from ( pT|) , 



T« = — [3^.(/Cl, ^2, A;i2)3^s(fc3, ^4, ki2) + (k2 ^ kg) + (k2 ^ k4)] 



(32) 



while the rest can be read off from (^) and (^3|). We emphasize that there are no other 
singularities in T: one can check that the logarithmic divergences appearing in intermediate 
formulas cancel out when one takes into account all contributions. 



Following Ref. ITQ] , we are going to compare our trispectrum to the trispectra of the local 
ms. The latter are obtained from the Ansatz in real space |2y, |2l[ 

C(X) = C. + \fNL{Cl - iO) + ^9NL{e, - 3(C,^)C,) , (33) 



*The product J^^J^r is regular in the limit ki2 0. 
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where (g is the Gaussian field and /nl and qnl are constants. One has |10 



G'4,ioc(ki, k2, ka, — — {Inl'Tioci + 9nlT\oc2) , 

where the two local shapes are 

= ^t*!- (35) 



1=1 



Let us quantify the strength of the non-Gaussianity. The standard estimator t^vL is related 
to the four-point function in the regular tetrahedron limit, ki = ki2 = = k, 



■■ki2=kid=k 



(27rfPj<5(^|^k,jlt^^. (36) 



Note that the standard definition involves in the right hand side. This is appropriate for 



the local Ansatz (p3D, and the sizes in the local models are 

^iVL = 2-16/^^, t^^f = l.OSgNL- 
On the other hand, the quantity we can directly calculate in our model is 



where Ve = l/(47r^), see (|T2D , and the size of the non-Gaussianity in 6 in our model is 
obtained from (||) and (H), 

C = 2.87/1^ . (37) 
The adiabatic perturbation ( is proportional to 6, namely, 

C = r| , (38) 

where 6 is the homogeneous background value of the scalar field, and r < 1 is a dilution 
factor, which is independent of k for both curvaton and modulated decay mechanism of 
conversion of the ^-perturbations into adiabatic ones. Therefore, in our model the size of 
the non-Gaussianity of the adiabatic perturbations is 
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As discussed in Section |T], there are no model-independent constraints on 9 and /i, so the 
intrinsic non-Gaussianity we study in this paper may well dominate over conversion-related 
one and be detectable. 

Let us point out, however, that in the conformal rolling scenario of Ref. [Q, the adiabatic 
power spectrum is itself proportional to K^. The reason is that G is the phase field. With 
our normalization, its background value is bounded from above, see Eq. (p!^ , 

and without fine tuning |0| ^ -n jh. Equation (|38|) then gives 

SO that ^ 

tNL < 2.87^ . 

Therefore, the intrinsic non-Gaussianity can be sizeable only for small dilution factor r. 
If the ^-perturbations are converted into adiabatic ones by the curvaton mechanism, the 
conversion-related non-Gaussianity is large at small r, /tvl ~ r^^ ||l6l, so its contribution to 



the trispectrum is very roughly of the same order as that due to the intrinsic non-Gaussianity. 
The analysis of the detectability of the intrinsic non-Gaussianity in this case deserves further 
study. The same remark applies to the bulk of the modulated decay models, where /atl ~ 
as well Ijl^. There is an exception, however |[T^: if the width of decaying particles r(^^) is 



linear in 6 or has the form r(^) = (70 + 71^)^ (which is more plausible from particle physics 
prospective), then /jvl, Qnl are roughly of order 1 even for small r ~ 71/70, and hence large 
tNL- We conclude that in the conformal rolling scenario, the domination and detectability 
of the intrinsic non-Gaussianity is possible, but not at all generic. 

Let us now turn to the shapes. To compare them, we set in what follows 

Inl = qnl = 1 ; 



then the sizes of the trispectra in all models are similar. 

Let us consider various limits of the shape function T. We use the nomenclature of 
Refs. |]10|, The first three panels in Figs. |l], |, |5 and ^ show the complete trispectrum T, 



the contribution of the singular part Ts and the regular part % in our model, respectively. 
The fourth and fifth panels show the two local trispectra TIoci and 7roc2- Note that vertical 
scales in the last two panels are different from each other and from vertical scales in the first 
three panels. The ranges of arguments in these figures are limited due to various inequalities 
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obeyed by the momenta. In particular, 



n4 



kf + kf- 2kiki < fci4 < \ kf + ki + 2kik^ 



The hmits we present are: 
1. Equilateral limit, ki = k2 



ks = k^. We plot in Fig. the trispectra as functions of 



ki2/ki and ku/ki. Clearly seen is the singularity T oc k^2^ k^^ in our trispectrum and in its 
singular part, as well as stronger singularity TToci oc k^2^ k^l in the first local trispectrum. 
As we pointed out in Section [1|, inflationary models produce trispectra without singularities 
at ki2 — )■ and/or ki^ — )■ [0, 0, ^ |T^. Thus, the singularity, which is due to the 
infrared enhancement of the modes 6p, is a distinctive feature of the conformal models. 
2. Specialized planar limit, ki = k^ = ki^ and 



k 



12 



kl + 



^2^4 
2^ 



k2k^+ J{Akl-kj){Akf 



1/2 



The trispectra are shown in Fig. |^ as functions of k2/ki and k^/ki. The structures along the 
diagonal are again due to the singularity, now at ki^ — )■ 0, which corresponds to k2 — ?■ k^. 
Note that our total trispectrum vanishes at the boundaries k2 = and /C4 = (this can be 



established analytically). The latter feature is similar to many inflationary models [|T0], |T5 
while it is absent for TToci and 7k)c2- 

3. Near the double-squeezed limit, ^3 = ^4 = k^. We show in Fig. ^ the combinations 

r 



i=l 



and in Fig. | the trispectra T themselves as functions of k2/ki and k^/ki. Clearly, our 
trispectrum is quite different from local ones in this limit. 



4 Conclusions 

By comparing the upper panels in Figs. |l], |], H, ||one observes that the most notable features 
of the trispectrum in conformal models are well captured by its singular part, which has the 
factorized form (|32|) . It is also clear that the trispectrum is substantially different from the 
trispectra of local forms. Furthermore, the comparison of our trispectrum with the trispectra 



given, e.g., in Refs. [0, ^ shows that our trispectrum is considerably different from the 
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trispectra inherent in inflationary models. Hence, the shape of the non-Gaussianity, together 
with the statistical anisotropy, is an interesting signature of the conformal mechanisms. 

As we observed in Section \i.2[ models employing conformal invariance for generating the 
flat scalar power spectrum are indistinguishable at the leading non-linear order. It remains 
to be understood whether one can discriminate between concrete models from this class, 
even in principle. In any case, it would be extremely interesting to learn (or rule out) that, 
in a certain sense, our Universe started out conformal. 

Note added. After this work has been published in arXive, the paper p4| appeared, 
where the results similar to those presented in Section ^ were obtained in much more general 
context. We are indebted to K. Hinterbichler and J. Khoury for making their paper available 
to us. 
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Figure 1: Complete trispectrum T (upper left panel), its singular part Ts (upper right 
panel), its regular part %■ (middle panel), trispectrum of local form TIoci (lower left panel) 
and trispectrum of another local form T\oc2 (lower right panel) in equilateral limit. See 
Section p.3| for definitions. Note that vertical scales in the lower panels are different, and 
neither coincides with the vertical scales in the three upper panels. 
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4 

Figure 3: Near the double- squeezed limit. Shown are the combinations T/(n ^i) the 

i=l 

same trispectra as in Fig. 
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Figure 4: Near the double- squeezed limit. Same as in Fig. ^, but for the trispectra T 
themselves. 
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Figure 5: s-channel diagram that contributes to the trispectrum. Dashed and sohd hues 
correspond to the two-point functions of the fields 6 and 5p, respectively. 

Appendix 

In this Appendix we perform the complete calculation of the tripsectrum. Notably, this can 
be done analytically. 

The computation is conveniently performed in terms of symmetric polynomials. For the 
s-channel diagram of Fig. ^ these are 



Zk„k, = K,{ml - I2K2 - kl^) , Xk„k, = kl^ + 2kUKl - 2K2) + UKIK^ - 3K^ . 



Ki = ki + k2 , K2 = kik2 , Pi = ^3 + ^4 , P2 = k^kA^ . 



In what follows, we encounter the combinations 




Q = [Kl-2K2-kl^-K2{Kl-kl^)e]sm{K,0 
- K^{Kf - 2K2 - kl2)icos{Kri) . 



The terms (|2^) give the following contribution to the four-point function: 
Gi Tf = ^ k2, ks, ki, ku) + 2 permutations 




i=l 



where, in self-explaining notations. 



00 
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and 



U{ki,k2,k3,ki,ki2;0 = j -^^(^J5/2{ki20y5/2{ki2X) - y5/2{ki20J5/2{ki2x) 





,fc4,fci2 (x) + Qk3 (0-Rfcl,fc2,fcl2(x) 

A tedious but straightforward calculation of the latter integral gives 

U = Uo + Ui, 

where 

uo = {2kUcos{K,om2ki,e - ml - ^k2 + ki,)) 

- 2Kiki2 sin(KiO(2A;t2e' + ^kl^ + 36/^2 - 9^^)^ + {k^, k2 ^ k^, k^}, (39) 

Ur = ^''''f^jP^''''''' [[{kW - 3) Mk,20 + Sk,2^cos{k,20] 

X [Ci((fCi - k^2)0 + Ci((i^i + A;i2)0] 
+ [{kU^ - 3) cos(A;i20 - 3A;i2e sin(A;i20] 

X [Si{{K, - kuK) - Si((i^i + ^12)0]) + {^1, k2 ^ ks, k,} . 

Note that Uq contains only trigonometric functions and powers of ^. So, the integration of 
Uq over ^ is cumbersome but straightforward^. On the other hand, Uj involves cosine and 
sine integrals Ci and Si, so the integration of it is tricky. We perform the latter integration 
by making use of the integral representations of the functions Si and Ci and changing the 



^At ^ ^ the original integrals converge. Nevertheless, some particular terms may produce divergences 
(note, e.g., that Uo contains cos(i^i^)/^). To regularize these divergences we integrate over ^ from a > 
and take the limit a — > in the end of the calculation. 
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order of integration, i.e., we first integrate over ^ and then integrate over 

oo 

32 f 



32 f cos{Kix)dx f ^^Qk3M,ki2iO^KuK2 



X 



X 



[(3 - kl2^'^) sin(/ci20 - ^ki2^ cosikuO] cos{ki2X) 
- [(3 - fciaC^) cos(A;i2C) + Sk^i sm{ki2i)] sin(/ci2x) ) + {ki, /ca ^ /cs, /i;4} 



Tfie inner integration over ^ is lengtliy but again straiglitforward. It again produces two 
types of contributions: the first one does not contain sine and cosine integrals. 
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16Xk,.k., f cos{Ka)dx I 2p^(^2^ _ 3P,)^eos(Pa) + + 12^2 - 3^^)^ cos(fc,2X) 



X 



+ [3Pf - 6P2 - k^S - 2P2X')] sin(Pix) j + {ki, k2 ^ ks, k^ , 
while the second one contains Si and Ci in the integrand, 
8Xki 



(40) 



U 



II 



Ci((Pi + Mx)-Ci((Pi-Mx) 



- log 



Pi + fcl2 

Pi - ki2 



+ sin(/ci2x) [Si((Pi + ki2)x) + Si((Pi - A;i2)x)] ) + {^1, ^2 ^ k^, k^ . 

(41) 



Nevertheless, both of these integrals can be evaluated analytically, the relevant formulas 
being 

00 

o / , X ,. f , _f^sin(aa;)Si(6x) 

SAa, b, c) = hm / dxe — ^ — L^—L + |5 ^ c| 



[e{b - a)Ll{b/a) + e(a - b)L2{b/a) - 2L2{-b/a) + L2(b/a)] + {6 ^ c} , 
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CAa, b, c) = lim / dxe~ 



. cos(ax)[Ci(6x) — Ci(cx) — log(6/c)] 



X 



-rUoK ^] lo 



[la)' 
be 



where L2 is the dilogarithm function. In this way we obtain 



U 



1 



hiM-^k-iM log 



2 \2 
12) 



{Pf - k 



{Kf - k 



2 \2 
12) 



V {K, + Pi 



+ 



j^^^f^,[Q2 + kl,g, + kt,g, 
1 



+ A;i2, Pi + A;i2, A - ^12) + - A;i2, Pi + A;i2, Pi - ^12) 

+ S,{K^ + A;i2, Pi + A;i2, Pi - ^12) - Si{Ki - ki^, Pi + ^12, Pi - k^) + {ki, k^ O k^, k^} 
The first two hnes here come from the integration of Uq, Eq. (|39|), and from (ffO), while the 



rest is due to (0). The notations are 

g2{ki,k2,k,,k,) = -12{K^ + P^)\Kl-AK2){P^-AP2), 

g4{ki,k2,k3,k4) = ?,2K2P2{Kl + K^Pi + Pl) 

- im{Ki + Pi){Kl + 2KlPi - miPl - 3Pf ) 

- IQK2{K^ + Pi) (Pf + 2PlK^ - SPi^i" - ml) 

- A{Ki + Pi)\3K^ - 2KlPl + Wt) , 

G^{k^,k2M,h) = -32^2^2-16(^1 + Pi) [K2(i^l+2Pi)+P2(Pl+2i^i)] 

+ 4(i^i + Pi)' (3^2 ^ Pi + 3Pi2) . 

We add the part of (^) that corresponds to the s-channel diagram of Fig. ^ and is given by 
(PSp, and obtain the overall contribution of the s-channel diagram: 

/l2 TT 



5(^1,^2,^3, ^4, 



12 



4 QO 

n 

i=l 



y{ki, k2, ki2)yik3, ^4, A;i2) + J{ki, /c2, ki2)J{k3, ^4, ^12) 



+ U{ki, k2, ^3, ^4, 



12 j 



(42) 



Together with crossing terms this yields our final result 



G4(ki,k2,k3,k4) 



i=l 



T{ki, k2, k3, /c4, A;i2, ku) = S{ki, /c2, k^, k^, k 



12 



+S{ki, /c3, k2, ^4, ^13) + 'S'(A;i, ^4, ^3, h, ^m) 



(43) 
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